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The single-freeze-out model with parametrized hypersurface and flow geometry is employed to 
analyze the transverse-momentum spectra of hadrons produced in the Pb+Pb collisions at the 
collision energy of t/snn = 2.76 TeV at the CERN Large Hadron Collider (LHC). With the notable 
exception for protons and antiprotons, we find a very good agreement between the model results and 
the data for the measured hadron species. The additional analysis of the HBT radii of pions helps 
us to select, from several different types of freeze-out studied in this work, the most realistic form 
of the freeze-out hypersurface. We find that discrepancy ratio between the model and experiment 
for the proton/antiproton spectra depends on pr, dropping from 2 in the soft region to 1 around 
p T = 1.5 GeV. 
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I. INTRODUCTION 

Thermal approach [T|-[25j has become one of the cor- 
nerstones of our understanding of ultra- relativistic heavy- 
ion collisions, allowing to explain the data on hadronic 
abundances collected at the AGS [IM1, SPS [H-GJ], 
and RHIC energies [USl. 

Amended with the single-freeze-out scenario and a 
proper choice of the freeze-out geometry [22T - |2~4j ] , it also 
properly describes the RHIC transverse momentum spec- 
tra [H, [^. collective flow [H, HI HI) femtoscopic ob- 
servables |30l - [32j , or certain two-particle correlation vari- 
ables, such as the charge balance functions dHHH. One 
obtains successful fits economically, with just a few ther- 
modynamic parameters, such as the temperature, T, the 
baryon and strangeness chemical potentials, hb and us, 
or in the extended approach the quark fugacities, j s and 
7 g , and several parameters describing the geometry and 
flow. For that reason in many popular approaches the 
thermal (or statistical) approach is used to model the 
hadronization stage of the heavy-ion reaction. 

Much to our surprise, the first measurements at the 
LHC [35| confirmed the validity of the thermal model for 
all identified particle species, but for the protons! This 
proton puzzle is indeed perplexing, as protons (and an- 
tiprotons) are basic products of the reaction and as such 
they need to be described in any successful approach. In 
this paper we argue that the puzzle is in some sense even 
deeper, as not only the abundances of all the other mea- 
sured particles are correctly reproduced, but also their 
transverse-momentum spectra can be described without 
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difficulty Thus, it is solely the protons (and antiprotons), 
where the model fails to describe these basic observables. 

In the present study we shall apply simple parameter- 
izations of the freeze-out geometry, similarly as in our 
early work [22T - I241 ] . which are motivated by the dynami- 
cal studies involving relativistic hydrodynamics l3ll . I36L - 
|4?| . In particular, we will explore the Cracow [22j and 
the Blast- Wave models QJ, [H |H, H| . The advantage 
of such an approach is that one can focus entirely on the 
freeze-out and flow geometry, avoiding intricacies of the 
earlier dynamical evolution. 

Our basic result is that the identified transverse- 
momentum spectra of all up-to-now measured particles 
at the LHC, with the important exception of protons, 
can be properly fitted, describing the Pb+Pb data at 
y/snn = 2.76 TeV in the soft regime. We also test the re- 
liability of the fits with the HBT correlation data, which 
allows us to select the most realistic freeze-out conditions 
out of several different types of freeze-out studied in this 
work. 



II. METHOD 

As mentioned above, the concept of the single freeze- 
out [22I - I24I |5C| allows for the uniform calculations of 
numerous soft hadronic observables. The simulations 
presented in this work have been carried out with 
THERMINATOR [HHl]. Monte-Carlo method allows for a 
simple inclusion of the experimental cuts and, therefore, 
for more realistic verification of the model predictions. 

At the extreme energies studied at the LHC, we expect 
that the midrapidity region contains equal numbers of 
baryons and antibaryons, hence the values of the chem- 
ical potentials used in this work are set equal to zero 
(pB = Us = A* J = 0, where /ii is the chemical potential 
related to the isospin conservation). Hence, we are left 
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with temperature, T, as the only independent thermody- 
namic parameter. A natural expectation, based on the 
shape of the T — fis freeze-out curve 0, 0] and confirmed 
by the analysis of the hadronic abundances at the LHC 
[35l ] , is that the temperature of the chemical freeze-out at 
the LHC should be very similar to the value determined 
at RHIC. Therefore, we use T chcm = 165.6 MeV [52|. We 
note that this value is close to the value of the transi- 
tion temperature of the crossover found at /ib = in the 
lattice simulations of QCD 

Having fixed the values of thermodynamic parameters, 
we test different parameterizations for the shape of the 
freeze-out hypersurface and for the form of flow at freeze- 
out to achieve the best description of the transverse- 
momentum spectra of charged pions and kaons. We stress 
that in our method only these two hadron species are 
used to determine the geometric and expansion parame- 
ters, since they are experimentally determined with the 
best accuracy. We fit the Cracow and Blast- Wave models 
(by the least squares method) to the experiment. Then, 
the spectra of other hadrons are predictions of our ap- 
proach, allowing for its verification. 

In the case of the Cracow model, we fit the value of the 
proper time at freeze-out and the transverse size of the 
fire-cylinder. In the case of the generalized Blast- Wave 
model, we fit the proper time and the magnitude of the 
transverse flow for three different values of the slope of 
the freeze-out curve in Minkowski space, controlled by a 
parameter A. The three different choices of A correspond 
to three different physical scenarios: if A is positive the 
freeze-out starts at the center of the system, if A is neg- 
ative the freeze-out starts at the edges, and for A = 
the freeze-out happens at constant (longitudinal proper) 
time in the whole volume. 

With the procedure outlined above, we find a very good 
agreement between the model results and the experimen- 
tal data for all measured hadron species excluding pro- 
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FIG. 1: Projection of the freeze-out hypersurface used in 
the Cracow model on the r — r plane (r = \/t 2 — z 2 and 
r = ^/x 2 + j/ 2 ). The parameters Tzt and r max are the two 
geometric parameters of the model. 



tons. Moreover, with the hadron spectra alone, no pref- 
erence for any of the employed freeze-out models can be 
found. A further analysis of the pion HBT radii indi- 
cates, however, that the freeze-out hypersurface used in 
the generalized Blast- Wave model with A = —0.5 leads 
to the best agreement between the data and theory. This 
confirms an earlier observation done for RHIC 3Q] . This 
type of freeze-out is also consistent with the hydrody- 
namic picture where the freeze-out starts at the edges of 
the system and continues inwards. 

III. CRACOW MODEL 

A. Definition of freeze-out conditions 

The Cracow model assumes boost-invariant and cylin- 
drically symmetric conditions at freeze-out. The freeze- 
out hypersurface is defined by the condition [22[ 

t 2 — x 2 — y 2 — z 2 — r 2 f = const., (1) 

and the fluid four-velocity is proportional to the space- 
time position 

u"=7(l,v) = — = — (l (2) 

T3f T3f \ t J 

The index 3 in the symbol rsf reminds us that the freeze- 
out takes place on the three-dimensional hypersurface of 
constant proper time, see Fig. [TJ 

Equations (flj and $2$ imply that the parame- 
terizations of the freeze-out hypersurface and the 
four-velocity field are connected. The longitudinal 
proper time, r = \Jt 2 — z 2 , and the transverse distance, 
r = \J x 1 + y 2 , may be expressed in terms of the trans- 
verse rapidity rj±, 

t = r 3f cosh?7 ± , 

t = r 3f sinhr^x. (3) 

We also use the parameterizations i = TC0sh?7|| and 
2 = Tsinh77|i where 7711 = 1/2 ln(t — z)/[t — z) is the 
space-time rapidity, and also x — r cos <f> and y — r sin <j) 
where <f> is the azimuthal angle. 

The calculation of the volume element of the freeze- 
out hypersurface in the Cracow model shows that it is 
proportional to the four- velocity (as in the original Blast- 
Wave model by Siemens and Rasmussen [Hj]), 

dS^ = u M r| f cosh(?7jJsinh(?7_L) drj± dm d<f) 

= u M T3f r dr dm dip- (4) 

The formula ((H) is used in the Cooper-Frye formula, 

- 27r Pl _ J d^(x)p»f(u„p»), (5) 

to generate hadronic states on the freeze-out hypersur- 
face defined by the condition (TTJ). In Eq. ([5]), the func- 
tion f{u v p v ) is the distribution function which is taken 
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FIG. 2: (a) Transverse-momentum spectra of pions (cir- 
cles), kaons (triangles), a nd p rotons (squares) for the cen- 
trality class c = — 5% [5y| . The Cracow model results 
(solid lines) generated with THERMINAT0R 2 [f^| were ob- 
tained for: T = 165.6 MeV, ^ B = 0, r 3f = 9.0 fm, and 
?"max = 11.4 fm. The geometric parameters were fitted to the 
spectra of pions and kaons, only, (b) The model ratio of the 
transverse-momentum spectra of A°'s and Ag's (solid line), 
compared with the data (triangles) for the same centrality 
class c = — 5% 57]. The values of the thermodynamic and 
geometric parameters are the same in the two parts. 



in the form of the Bose-Einstein or Fermi-Dirac distri- 
bution for bosons and fermions, respectively, and p^ is 
the particle's four- momentum expressed by the rapid- 
ity y and the transverse-momentum p± (the quantity 
mj_ = m? + p 'j_ is the transverse mass), 



(mj_ coshy,p_i_ cos 4> p ,px sin0 p , m± sinh y) . (6) 



The generation of particles is performed with 
THERMINAT0R 2 [H3] which includes all known hadrons 
containing u,d and s quarks (THERMINATDR 2 has the 
same particle input basis as SHARE (Hf). THERMINAT0R 2 
simulates also decays of resonances, hence, the final spec- 
tra consist of primordial and secondary contributions 
(the primordial particles are emitted directly from the 
fireball, while the secondary particles come from the res- 
onance decays). 

The Cracow single- freeze-out model has altogether four 
parameters: two geometric and two thermodynamic ones. 
The two geometric parameters are rgf and r max , while 
the two thermodynamic parameters are temperature, T, 
and baryon chemical potential, /i b ■ The thermodynamic 
parameters define the local equilibrium distribution func- 
tions f(u v p v ) in Eq. (J3J) - 



Since the values of T and \ib follow from the analyses 
of hadron abundances, in the present work we may use 
the results of the previous studies where the conditions 
for chemical equilibrium at the LHC have been studied, 
58]. On more general grounds, in the central region of 
the heavy-ion collisions performed at the LHC energy, we 
expect equal numbers of baryons and antibaryons which 
implies /is = (the other chemical potential should also 
vanish). This leaves us with temperature as the only 
independent thermodynamic parameter. Expecting that 
the freeze-out temperature is the same as that found for 
RHIC HI we use the value, T = 165.6 MeV. 



B. Comparison with the LHC data 

In Fig. [5] (a) we show the experimental transverse- 
momentum spectra of pions (circles), kaons (trian- 
gles), and protons (squares) 56] for the centrality class 
c = — 5% 1 . The Cracow model results generated 
with THERMINAT0R 2 [H3] are denoted by the solid lines. 
The parameters used in the model calculations are: 
T = 165.6 MeV, \i B = 0, r 3f = 9.0 fm, and r max = 
11.4 fm. The values of the thermodynamic parameters 
have been treated as external parameters and the ge- 
ometric parameters were fitted (with the least squares 
method) to the spectra of pions and kaons, only. 

In Fig. [3] (b) we show the model ratio of the transverse- 
momentum spectra of A°'s and ^5's (solid line), com- 
pared with the data (triangles) for the same centrality 
class c = — 5% [13] • The values of the thermodynamic 
and geometric parameters are the same in the two parts 

of Fig. m 

We observe very good agreement between the data and 
the model results for pions, kaons, and A°'s in the soft 
momentum region, p± < 3 GeV. On the other hand, the 
model results for protons overpredict the data. This is 
an expected result, since already the results for the abun- 
dances have indicated that the thermal models used in 
the grand-canonical version cannot predict correctly the 
ratio of pion and proton abundances [35[. 

In Fig. [3] (a) we show again the transverse-momentum 
spectra of pions (circles), kaons (triangles), and protons 
(squares) but now for the centrality class c = 10 — 20% 
[5 61 ] . The geometric parameters used in the model calcu- 
lation are: T3f = 7.4 fm, and r max = 9.6 fm. Similarly to 
the case of central collisions, the geometric parameters 
were fitted to the spectra of pions and kaons, only. In 
the natural way, they are smaller than those found in the 
case of the most central collisions. 



1 All the data and the model fits refer to Pb+Pb collisions at 
y / %J = 2.76 TeV. Moreover, to be consistent with the exper- 
imental procedure |57( . the spectra of protons are corrected for 
the decays of A°'s and TP's. Similarly, A°'s are corrected for the 
decays of S — 's. 
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FIG. 3: (a) Transverse-momentum spectra of pions (circles), 
kaons (triangles), and protons (squares) for the centrality 
class c = 10 — 20% [5q|. The Cracow model parameters used 
in the calculation: T = 165.6 MeV, \ib = 0, T$t = 7.4 fm, 



and 



= 9.6 fm. Again, the geometric parameters were 



fitted to the spectra of pions and kaons, only, (b) The model 
transverse-momentum spectra of A°'s, H~'s, and f2~'s, com- 
pared with the experimental results for H~'s and fi _ 's for the 
centrality class c = — 20% [3^|. The values of the thermo- 
dynamic and geometric parameters are the same in the two 
parts. Since the geometric parameters were fitted to the cen- 
trality class c = 10 — 20%, the model results for hyperons are 
slightly above the experimental data. 



In Fig. 131(b) we show the model transverse-momentum 
spectra of A 0, s, S~'s, and _ 's. They are compared 
with the experimental results for S _ 's (squares) and f2 _ 's 
(triangles) for the centrality class c = — 20% 57]. The 
values of the thermodynamic and geometric parameters 
are the same in the parts (a) and (b). Since the geometric 
parameters were fitted to the centrality class c = 10 — 
20%, the model results for hyperons are slightly above the 
experimental data but their p± dependence is consistent 
with the data. 

Similarly to the central collisions, we see again that the 
agreement between the data and the model predictions 
is very good except for the protons. 

In Fig. [4] (a) we show once again the transverse- 
momentum spectra of pions (circles), kaons (triangles), 
and protons (squares) but this time for the centrality 
class c = 30 — 40% (56[. The Cracow model param- 
eters used in the calculation are: — 5.9 fm, and 
r mal = 7.25 fm. They were again fitted to the spec- 
tra of pions and kaons, only. The decreasing trend of raf 
and r max with increasing centrality reflects the decrease 
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FIG. 4: (a) Transverse-momentum spectra of pions (circles), 
kaons (triangles), and protons (squares) for the centrality 
class c = 30 — 40% The Cracow model parameters used 
in the calculation: T = 165.6 MeV, /is = 0, Tat = 5.9 fm, and 
rmax = 7.25 fm. Once again, the geometric parameters were 
fitted to the spectra of pions and kaons, only, (b) The model 
ratio of the transverse- momentum spectra of A°'s and Kg's, 
compared to the data for the centrality class c = 20 — 40% 
[57j . (c) The model transverse-momentum spectra of A°'s, 
H~'s, and n~'s, compared with the experimental results for 
H" 's and fT 's for the centrality class c = 20 - 40% The 
values of the thermodynamic and geometric parameters are 
the same in the three parts. 



of multiplicity 

Figure 2] (b) shows the model ratio of the transverse- 
momentum spectra of A°'s and Kg's, compared to the 
data (triangles) for the centrality class c = 20 — 40% 
[I?}, whereas Fig. 2] (b) shows the model transverse- 
momentum spectra of A°'s, S~'s, and Q~'s, compared 
with the experimental results for S~ 's (squares) and fi - 's 
(triangles) for the centrality class c = 20 — 40% [35] . The 
values of the thermodynamic and geometric parameters 
are the same in the three parts. Similarly to Fig. [31 the 
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model hyperon spectra are slightly above the data, which 

may be explained by different centrality classes analyzed ^ c [%] r [ J 

in the presented comparison. 

Our comparison of the transverse-momentum spectra 
obtained within the Cracow model with the experimental 

data may be summarized with the statement that the " 

model describes well the spectra of all measured hadrons 10-20 8.3 0.375 

except for protons, in the p± range up to 3 GeV. 0.5 30-40 6.5 0.425 



IV. GENERALIZED BLAST- WAVE MODEL 

A. Definition of freeze-out conditions 

Probably, the most popular parameterization of the 
freeze-out hypersurface is the Blast- Wave model [3, [H[ ■ 
In its standard form, the model is boost-invariant and 
cylindrically symmetric (similarly as the Cracow model 
discussed in the previous Section). The Blast- Wave 
model uses the assumption that the freeze-out happens 
at a constant value of the longitudinal proper time 



T~2f — const. 



(7) 



In order to get a broader applicability, we generalize this 
condition to the formula 



t = r 2 f + Ar, 



(8) 



where r 2 f and A are constants, and A describes the slope 
of the freeze-out curve in the Minkowski space, see Fig. [5] 
With A > (A < 0) we may consider the freeze-out 
scenarios where the outer parts of the system freeze-out 
later (earlier). Of course, with A = we reproduce the 
standard Blast- Wave parametrization (J7J)- 

In the generalized Blast- Wave model we find compact 
expressions for the argument of the equilibrium distribu- 



T 2 f 



freeze-out curve 
t = t 2 i + A r 



FIG. 5: The freeze-out curves considered in the modified 
Blast- Wave model H. 



0.0 0-5 9.9 0.45 
0.0 10-20 8.2 0.43 
0.0 30-40 6.2 0.46 



-0.5 0-5 
-0.5 10-20 
-0.5 30-40 



10.4 0.46 
8.9 0.475 
6.9 0.58 



TABLE I: Optimum choices of the parameters r max and vt 
for three fixed values of the parameter A and for different 
experimental centrality classes. The fits have been performed 
to the transverse-momentum spectra of pions and kaons only. 
The fit has been constrained by the condition T2f/r m ax = 1. 



tions functions, 

[mj_ cosh(?7|| —y) — vj_(r)p^ cos(</> - tf> p )] 



u v -p 



and for the Cooper-Frye integration measure, 

= (r 2 f + Ar)r [m_L cosh(f?|| - y) 
-Ap± cos((/> - <0 . 



(9) 



(10) 



Equations © and (fl~U)) should be used in the Cooper-Frye 
formula ([5]). 

The user of THERMINAT0R 2 [Hf may choose different 
r-profiles of the transverse flow v±(r). In this work we 
use the following option 



r/r 



v T + r/r n 



(11) 



where vt is the parameter controlling the strength of the 
transverse flow. In this version of the Blast- Wave model, 
we have four parameters: A, T2{, r max , and vt- 



B. Comparison with the LHC data 

The freeze-out conditions defined by the generalized 
Blast- Wave model with different values of A were studied 
in Ref . 30] . One of the conclusions of this work for RHIC 
is that the transverse-momentum spectra of hadrons can 
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FIG. 6: Same as Fig. [2] but the data are compared to the 
results obtain with the Blast-Wave model with A — —0.5. 
The data are taken from [56l. [57j. 
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FIG. 7: Same as Fig. [3] but the data are compared to the 
results obtain with the Blast-Wave model with A — —0.5. 
The data are taken from [35l. |56^. 

be well reproduced for a wide range of the parameters A, 
however, the HBT radii are quite sensitive to the space- 
time profile of the freeze-out hypersurface, i.e., to the 
specific choice of the parameter A. 
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FIG. 8: Same as Fig. [4] but the data are compared to the 
results obtain with the Blast- Wave model with A — —0.5. 
The data are taken from [35], IBlil. [57l]. 

Inspired by Ref . [3(| , we have performed the analysis of 
the LHC data in the similar way. At first, we have chosen 
three values of A (A — 0.5, 0, —0.5) and for each of these 
values we have found the optimal geometric parameters. 
In order to reduce the number of independent parameters 
we have fixed the ratio t%{ /r max to unity. In the next step, 
for each value of A we used the optimal choice of 7' max 
and vt to calculate the HBT radii. 

The optimal choices of the parameters r max and vt for 
three fixed values of the parameter A and for different 
experimental centrality classes are given in Table [U We 
have found, as suggested by Ref. [3(| , that the transverse- 
momentum spectra are equally well described for differ- 
ent choices of A, if the other parameters are properly 
chosen. 

In Figs. [BHH] we show our fits done with A = —0.5. 
The three figures correspond to Figs. [2HH presented ear- 
lier in the context of the Cracow model. We observe 
again rather good agreement between the model predic- 
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FIG. 9: The pion HBT radii for the most central collisions 
obtained in the Cracow model (solid lines) and compared to 
the LHC data from [59(]. The parameters are the same as in 
Fig.H 



FIG. 10: The pion HBT radii for the most central collisions 
obtained in the generalized Blast- Wave model with A — —0.5 
(solid lines) and compared to the LHC data from [59]. The 
parameters are the same as in Fig. [6] 



tions and the data. Noticeable discrepancies can be ob- 
served in the ratio of the transverse-momentum spectra 
of A°'s and K^'s for the centrality class c = 20 - 40%, 
see Fig. |S] (b). However, these differences are the largest 
(about 20%) for the momenta reaching 3 GeV, i.e., in the 
region where we expect the thermal approach to break 
down. Moreover, since the thermal approach uses many 
simplifying assumptions, the agreement within 20% is 
usually regarded as quite satisfactory. 

We do not show here our results obtained for the cases 
A = 0.5 and A = 0. They are very similar to those 
presented in Figs. [2HH an d HHH We note that the shape 
of the freeze-out hypersurface in the Cracow model is 
similar to that used in the generalized Blast- Wave model 
with A = 0.5 (along the freeze-out curve the proper time 
grows with the distance from the center). This suggests 
that the results obtained with the two models should be 
similar and the values of the optimal close to each other. 
Indeed, the proper time used in the Blast- Wave model, 
T2f , is to a good approximation an average of the proper 
time and the transverse size used in the Cracow model, 

7~2f ~ (^*3f ~t~ ''max )/2- 



V. HBT RADII 

Our results presented in the previous Sections indicate 
that, except for the protons, the transverse-momentum 
spectra of different hadronic species may be well re- 
produced in the thermal approach. Moreover, different 
freeze-out conditions may lead to very similar spectra, as 
it has been demonstrated by our results obtained with 
the two different versions of the thermal model. In this 
Section we present calculations of the HBT radii done 
in the Cracow model and the Blast- Wave model with 
A = —0.5. It turns out that the use of an additional 
observable may, in the considered cases, select the most 
appropriate version of the freeze-out model. 

Figure [9] shows the pion HBT radii for the most cen- 
tral Pb+Pb collisions, c = — 5%. Theoretical results 
obtained in the Cracow model (solid lines) are compared 
with the LHC data taken from Ref. [5{| (i? ut _ circles 
in part (a), i? s ide _ squares in part (b), and i?i on g ~ tri- 
angles in part (c)). The radii are presented as functions 
of the transverse-momentum of the pion pair. They have 
been calculated in THERMINAT0R 2 with the help of the 
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two-particle method (without Coulomb corrections). The 
parameters used in the simulations are the same as those 
used to obtain the spectra shown in Fig. [2j We find that 
i? s idc is well described, i?i on g is a bit to large, while i? ou t 
is clearly underprcdictcd. 

In Fig.[lj)]we show analogous results for the Blast- Wave 
model with A — —0.5. With much improved agreement 
of i? ut with the data, the general consistency between 
the model calculations and the data has been significantly 
improved. 

VI. CONCLUSIONS 

The single-freeze-out model with parametrized freeze- 
out hypersurfaces and flow has been used to ana- 
lyze the transverse-momentum spectra of hadrons pro- 
duced in Pb+Pb collisions at the collision energy of 
^/snn = 2.76 TeV at the LHC. Except for protons, we 
find a proper agreement between the model results and 
the data for all measured hadron species. The additional 
analysis of the HBT radii of pions suggests that the real- 
istic freeze-out conditions should correspond to an earlier 
freeze-out of the edges of the system, as suggested by a 
typical form of the hydrodynamic expansion. 

Thus, the LHC proton puzzle is a certain sense deeper: 



it appears not only in the simple thermal approach 
where the abundances are calculated, but also in an 
extension where the pr-spectra can be computed. Of 
course, with the mismatch on abundances, which is a 
PT-integrated measure, we need to find a mismatch in 
the proton/antiproton spectra. We note, however, from 
Figs. [2] and [6] that the model proton/antiproton spec- 
tra are above the data only in the soft region below 
Pt ~ 1.5 GeV, while the harder part is in agreement. 
Therefore the proton puzzle is clearly related to the soft 
physics. We also note that the ratio of the model to 
experiment is pT-dependent, dropping from a value of 
about 2 at the 1ow-pt values to 1 at px ~ 1.5 GeV. 
Thus a simple rescaling of the proton/antiproton spec- 
tra, if found in some treatment, would not do the job. 
This issue, crucial for the thermal approach to relativis- 
tic heavy-ion collisions, requires a further study. 
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